Superbosonization formula and its application to random matrix theory 
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Starting from Gaussian random matrix models we derive a new supermatrix field theory model. 
In contrast to the conventional non-linear sigma models, the new model is applicable for any range 
of correlations of the elements of the random matrices. We clarify the domain of integration for the 
supermatrices, and give a demonstration of how the model works by calculating the density of states 
for an ensemble of almost diagonal matrices. It is also shown how one can reduce the supermatrix 
model to the conventional sigma model. 
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I. INTRODUCTION 



' The non- linear supersymmetric sigma model (NLcrM)! has been very fruitful in describing quantum systems with 
c/2 , quenched disorder—. Derived originally^ from a model of noninteracting electrons in a Gaussian random potential, 
^ ■ NLcrM was later constructed in its zero-dimensional version directly"^ from the Wigner-Dyson ensembles'* of random 
^ , matrix theory. This gave the first demonstration of equivalence between these two problems in the low-frequency 

■ domain. Furthermore, even for problems where the conventional Wigner-Dyson random matrix theory is sufficient 
^ _ (see, e.g.— 1^1^) NLctM turned out to be a competitive method of calculation compared to the standard approach based 
^ . on orthogonal polynomials-. Let us mention parametric level statistics^ as a beautiful example of new results which 

were obtained in the field of classical RMT using NLcrM. 

The most significant applications of NLcrM, however, are in the realm beyond classical RMT, such as corrections 
Q to the Wigner-Dyson theory in systems with diffusive dynamic o^'^'^^ , and random matrix ensembles with a proba- 
cy , bility distribution not invariant under the full group of unitary transformations. The domain of application of such 
random matrix ensembles is rather wide. For example, band random matrices are equivalent to thick disordered 
I ' wires^iii, where strong localization was first established using NLcrM. Another example is the Critical Power-Law 
. Band Random Matrix (CPLBRM) ensemble^^. An important property of eigenvector statistics in this ensemble is 

■ multifractalit y^^'^'^i^^ , which makes CPLBRM a good model to qualitatively (and in many cases quantitatively) de- 
scribe the critical states near the Anderson transitio n^'^i^^ . In the analysis of all these models NLcrM was an extremely 

On useful guide, although limited as NLcrM is not always exactly equivalent to the model being studied. The reason for 
, this limitation is the saddle-point approximation which is used in deriving NLcrM. In the case of random band matrices 
t — ■ the band width B has to be large in order for the saddle-point approximation to be valid. For CPLBRM this implies 
' weak multifractality, while for random band matrices it means that the localization length is much larger than the 
lattice scale. In both cases localization effects are in a certain sense weak. 

There are several important problems which lie beyond the applicability of NLcrM. These include the properties 
of one-dimensional disordered chains and the problem of critical eigenfunctions and level statistics in systems with 
strong multifractality. The Hamiltonian of the first example is just a tridiagonal random matrix, while a representative 
^ , model for the second example is CPLBRM with a small bandwidth i? ~ 1 . Both of these are examples of a class of 
' Gaussian ensembles of almost diagonal random matrices which were recently defined and studied in RefsJ^'i^iii^'i^. In 
order to attack these problems by the supersymmetry method, one needs an exact representation of the corresponding 
correlation functions in terms of a supermatrix field theory. 

Such a representation called 'superbosonization' has recently been suggested in Refsi^i^i. Similar ideas had been 
put forward earlier by Lehmann, Saher, Sokolov, and Sommersj^^ and by Hackenbroich and Weidenmiilleri^. Here 
we give a decent derivation of the key formula of super-bosonization, including such aspects as integration measure 
and domain of integration. We also work out a simple application of the formula involving Gaussian random matrices. 

Given any function f{tp®^) of the tensor product of a super-vector tl) and its conjugate ^, the super-bosonization 
formula essentially converts the integral of / over to an integral over a supermatrix Q akin to ® {p. Thus it 
provides a direct way to rewrite in terms of Q the so-called ^/^-functional which arises after ensemble averaging the 
initial free-field representation of disordered onc-particle systems. We wish to stress that the conversion from ip to Q 
is exact. It is also non-trivial, as the number of independent variables in the supermatrix Q may be smaller or larger 
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than the number of independent variables in the initial super-vector representation. We note that a workabl o-'^^i^ 
integration formula which converts integrals over tp into integrals over some supermatrix Q (different from Q) can 
be constructed by a mere change of variables^^ without changing the number of variables. What we describe in this 
paper is not simply a change of variables. 

For ease of discussion, we consider disordered one-particle Hamiltonians in the form of a random Hcrmitian matrix 
H with independent, Gaussian-distributed entries Hij = 1, . . . ,N) characterized by a certain variance matrix 
Cij = c{\i — j\) = The special case of classical Wigner-Dyson random matrix theory is obtained by setting 

~ j|) = const. While this has the largest possible symmetry group, U(-/V), our ensembles with non-constant 
c(|i — j|) generically have symmetry group U(l)^. From the perspective of possible applications, the main outcome of 
this paper is an exact field-theoretical representation of correlation functions for the generic case of U(l)^ symmetry. 

To exert better control over the mathematics, and also with a view to certain applications, we consider random 
matrix models with a local gauge symmetry U(n)^/". The random matrix Hamiltonians with matrix elements Hf^ 
in such models have a block structure with z,j — 1, ...,M := N/n and a,b = 1, ...,n. The probability distribution 
function is of the form 

(AI n 
ij — l a,b—l 

Setting M = 1, n — N gives the classical Gaussian Unitary Ensemble (GUE) of Wigner and Dyson, while the case 
n = 1, M = N corresponds to a generic Gaussian random matrix model. 

It turns out that the precise form of the superbosonization formula depends in a crucial way on how n compares 
with the number of commuting components q in the supervector ip. The formula agrees with the naive expectations^!^ 
a q < n (as is the case, e.g., for the mean density of states (DoS) of a generic Gaussian Hermitian random matrix 
ensemble, where q ^ n = 1) but has a different form ii q > n (which is of relevance to the DoS two-point correlation 
function of a generic Gaussian Hermitian random matrix ensemble, where q — 2, n — 1). Establishing and interpreting 
this peculiar fact is the main message of this paper. 

The structure of the paper is as follows. In Section[lllwe derive the super-bosonization formula and the corresponding 
supermatrix field theory. Our detailed treatment makes it clear that the formula changes as one goes from one-point 
functions (i.e., the mean density of states, DoS), to two-point functions (e.g., the DoS-DoS correlation functions), 
to three-point functions, and so on. In Section lllll we use the formalism to compute corrections to the DoS of the 
Gaussian ensemble of almost diagonal random matrices. These were obtained earlier in Ref by another method. In 
Section HVl we show how the new supermatrix field theory can be reduced to the conventional NLcrM. In Section fVl 
we give a summary of the main results and outline some open problems. In particular we discuss the possibility of a 
unified description of the cases q > n and q < n . 

II. THE SUPERBOSONIZATION FORMULA 

The first few steps of the method presented here are similar to previous derivations of the non- linear sigma model for 
random matrices^i^. However, while the usual scheme employs a Hubbard-Stratonovich transformation after the initial 
step of averaging over the random matrix, here we will follow a different route avoiding the Hubbard-Stratonovich 
transformation. 



A. Review of basic steps 

We begin with a generating function Zq [J] , 

Zo [J] = J D^PD^P exp (-i E^ ^ sUij ^j) , % = H,j - 5,, {£ + .h) , (1) 

where ipi are supervectors (with components arranged as column vectors), 'ipi are the conjugate row supervectors, Hij 
are the matrix elements of our random matrix iJ, Ji is the source field, and £ is a diagonal matrix containing the 
energy parameters of the problem. The diagonal matrix s with entries ±1 keeps track of whether the Green's functions 
to be computed from the generating function ([1]) are retarded (i.e., have an energy parameter with positive imaginary 
part) or advanced (negative imaginary part). A scalar product or sum over components 'ipiSipj = 'Ylix'^f^^'^^ 
in the exponential is understood. The number of components of the supervectors ij^i is chosen according to the 
physical quantity under consideration. In the most general situation each ipi has p anti-commuting and q commuting 
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components. To calculate such quantities as the two-level correlation function one sets p ^ q = 2 . In the case of the 
density of states it suffices to use p — q — 1 . 

In order to make the presentation self-contained, we now give more details concerning the supervectors ^pi and the 
matrix s . The definitions are the same for every site, so we temporarily drop the site index i from our notation. We 
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n have — ygj ■0 = (x where 

, X - (x'* ■ ■ • x^*) , s = {s^* . . . s^*) , 

and and S"^ are anti-commuting and commuting variables, respectively. S"^* is the complex conjugate of the 
complex variable S'^ . The integral ^ is carried out w.r.t. the Berezin superintegral form 

D^D^ = (27r)-P dSdSd^d^ , dSdS = [] 2 dmciS^ d3m{Sn , d^d^ = J] ^ ^ , 

u=i tj.=i ^ ^ 

with the domain of integration for each commuting variable S'^ being the complex numbers, C . The integration 
measure for the multi-site problem ^ is of course D-ipDip = D^JiDf/ji (index i reinstated). To specify the sign 
matrix s , let us agree that we want <?+ retarded and q- advanced Green's functions (where g+ + g_ = q)- Then we 
may choose s to be the diagonal matrix 

s = diag(lp,l5+ ,-19-) ■ 

(Note that there is absolutely no gain from introducing unnecessary minus signs in the fermion-fermion sector.) 
The next step is to take the average over the random matrix using the Gaussian distribution 



)h^M j dH{...) exp (-i , I' l^'^o) : 



where the coefficients are real, non-negative, and symmetric under exchanging i ^ j , and M is the normalization 
constant ensuring {1) h — ^ ■ Of course H*^ = Hji by Hermiticity. The integration measure is 

dH := [| dHkk n 2 dmz{H,j) d3m{H,,) , 

k i<j 

where we are using a normalization convention which will be in force throughout the paper. 

An equivalent way of characterizing such a Gaussian distribution is by means of its Fourier transform with respect 
to some arbitrary commuting parameters Kij : 

exp (^-i _ HijKj^^ ^ = exp ^ . ^ C^jKi^Kj^ . (2) 

The supervector integral representation ([T]) of the generating function Zi [J] allows us to average immediately over 
the random matrix H . Using Eq. ([2]) for Kji = tpi s -ipj , the average Z[J] of the generating function is 



Z[J] EE {Zo[J])h = J D^D^l^ exp (i i^^s{E + J{)^, " ^ I^. , ^^^^ V^jOlV^j^ V'^) 



(3) 



Up to now, everything was more or less standard. In the next step, however, we shall depart from the standard 
method^, where one proceeds by making a Hubbard-Stratonovich transformation. Instead, we will carry out a certain 
change of integration variables. The change-of-variables formula to be introduced below bears some resemblance to 
the well-known Weyl integration formula, by which one reduces the integral over a conjugation-invariant function of 
a matrix to an integral over the eigenvalues of the matrix. 

From the cyclic property of the supertrace, one has the identities 

STr (V'i ® '4>i)s ^ ^iS^pi , STr (ipi (8) -1^1)3 {ipj (g) ^j)s = {{pi s ipj){ipjS ipi) . 
Reading these identities backwards, one reorganizes the exponential of the integrand in ([3]) as 

i ^ . STr (V-i (g) Vi,)s (£ -t- J,) - i ^ . . C^j STr (g) 'iP,)s (ipj (g i^j)s , 
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which prompts a change of integration variables 

i^z<»i^z^Qt , (4) 

where Qi is a supermatrix. Note that while the sign matrix s does appear in the integrand, it will have no influence 
on the definition of the integration variable Qi . Thus the domain of integration for Qi will be independent of the 
types of Green's functions (retarded or advanced) which are being generated. This is a major difference from the 
traditional method using the Hubbard-Stratonovich transformation. 



B. Heuristic approach 

The change of variables tl^i ® il^i —> Qi works in the same way at every site i . Let us therefore focus on a fixed site 
and simplify the notation by dropping the index i for now. 

In past work on the subject the desired change of variables @ was brought about by an unprecedented and 
unexplained method, in the context of supermatrices, involving the (5-function. The logic went something like this: 
to transform a Berezin superintegral such as J DrpDip fiip (g) rp), you insert 

1 = j DQ 5{Q - ® i)) 

under the integral sign. Here Q is a supermatrix with an unspecified domain of integration. Then you reverse the 
order of integration of Q and ip to write 



Di^Di; /(V- ® V) = j DQ J{Q) f{Q) , 
J(Q) = / D^Di: 5{Q-ij;®^) = SDet(Q) 



This kind of manipulation with the (5-function requires justification, which has never been provided. To give a 
simple indication of what the issue is, consider the fermion-fermion sector for the case of the two- level correlation 
function, i.e., put p — 2 and q = , in which case we are dealing with 2x2 matrices 

(^xV* xV*\ Q_(Q'' Q'' 
^®^-UV* xV*j ' ^-\Q^' Q^' 

For concretcness let us say that /(■(/; ® iji) — f{x 'S) x) = X'^X^*x'^X^* ■ Now, in attempting to decide what function 
f{Q) is to be placed in the integrand we encounter an ambiguity. If we group the anti-commuting variables as 
/ = (x^X^*)(x^X^*) we may be inclined to choose f{Q) = Q^^Q^^ , but it is equally valid to reorder the anti- 
commuting variables as / = — (x^X^*)(x^X^*): which would suggest f{Q) = —Q^^Q^^. Perhaps a linear combination 
f{Q) = Q"Q22 _ ^^giiQ22 ^ q12q21>, ^^j^-j^ t G C) IS the good choicc to make. All of these do the necessary job of 
returning the given function / — x^x}*X^X^* on making the substitution Q ^ X® X- Thus, unless the matrix Q is 
constrained by Q^^Q^^ -|- Q^^Q^^ = , an arbitrary parameter t enters into the calculation. 

At this stage it must be recalled that in mathematics there exists no such thing as the S-function. What does exist 
and can be made sense of is the 6- distribution, say Sx , which is the linear functional that evaluates test functions at 
the point x : 



f^S4f] = I S{x - y) f{y) dy f{x) 



In our example, however, we are dealing with the element / — x}'x}*X^X^* of ^ Grassmann algebra - not with a 
function defined on points. It is unclear how we can evaluate a (5-distribution over these variables. In other words, 
the meaning behind 



m-x®x) = HQ'' - x\n HQ'' - rri m'' - x'x'*) m'' - x'x'*), 

viewed as a 'function' of the Grassmann variables, is undefined. No meaning exists in current mathematics. 

In the simple case under consideration, it is not difficult to work out everything by hand and achieve control of the 
situation. Let us take the integration domain for the 2x2 matrix Q to be the unitary group U(2) - or any other 
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4-dimensional closed submanifold of GL2(C) in the same homology class as U(2). Then if dvohiQ) denotes a Haar 
measure on U(2), it can be shown that 

/ dvoi2(Q)j(Q)(Q"g22 + gi2Q2i)^o, j(g) = Dct-i(Q) . 

So, although bosonization oi f — X^X^*X^X^* is an ambiguous procedure leading to any one of the one-parameter 
family of functions Ft{Q) = Q^^Q^^ — t{Q^^Q^^ + Q^^Q^^), the ambiguity, at least in this case, disappears at the level 
of integration provided we use J{Q) = Det~^{Q) and the integration domain U(2). 

This does not solve the problem of what to do when the Grassmann variables are more numerous or when considering 
the full supersymmetric situation which is rendered more complicated by the presence of commuting variables. It is 
clear that the same ambiguity occurs in the general case: for any given function /(■(/' ^ there exist many choices 
of function F{Q) such that F{Q)\q_^^^^ is equal to f{tp (g) -0). Given all these possible functions F{Q) with which 
the superbosonization formula 

J D^D^f{^ ®^)^ J DQ J{Q) F{Q) 

appears to be true, we may think we can use any one of them, or we may attempt to impose some constraint on 
Q. Fortunately, recent mathematical work2^ has given a complete solution to this problem for a restricted range of 
parameter values. In the next subsection, we will present a summary of these results. Afterwards, will give a detailed 
proof of the superbosonization formula for a special but important case. 



C. Rigorous result 

In order to control the mathematics, let us assume that there is a block or granular structure in the variances Cij . 
By this we mean if i is a multi-index i = (/. a) where a = 1, . . . , n , then 

Cij — Cl_a;I',a' — CH'Saa' i 

independent of a, a'. If so, then after grouping the terms appropriately, the integrand depends only on the sums 
^^^-^ ^pi^a i^i' ,a ■ In this way the integer n (sometimes referred to as the number of orbitals) is introduced as an 
additional parameter of our problem. Results stated in this subsection are valid only in the range n > q . Later, we 
will suggest how to recover the important case of n = 1 when q > 1 . 

To keep the notation simple, we return to using a single index i as a label for our supervectors ipi . In keeping with 
the discussion above, we consider functions / of the sum ''Pi®^i ■ Our goal is to transform (by superbosonization) 
the integral of such a function, 

j f = j D4'D4, f{^^^^® i^^j , D4,D^ = JJ Di>,Dil^, . (5) 

Let us address the general case of supervectors V^,; — that havep anti-commuting and q commuting components 

(for Xi and Si respectively). Then X^ILi ^« ® ^« corresponds to a supermatrix 



Q = 



A a 
T B 



where the blocks A and B are square matrices of size p x p and q x q with commuting variables as entries; while a 
and r are rectangular matrices of size p x q and q x p respectively and these have anti-commuting entries. 

Let F now be any function of the supermatrix Q so that on making the substitution Q ipi ''Pi the function 

F{Q) becomes equal to the given function f{J2 4'i®4'i)- We also assume that n> q. Then we claim that the following 
equality of integrals holds: 

J D^DiP / ff^ V'. ^ = ^oiu(l+^p-g) SDet"(Q) F{Q) , (6) 



provided / decreases sufficiently fast at infinity so that the integral on the left-hand side exists. 
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We now define all terms on the right-hand side of the identity ([6]). The domain of integration D is the unitary 
group U(p) for A and is the positive Hermitian q x q matrices, Herm^(<7), for B . Thus D — U(p) x Herm^(<7). The 
Berezin superintegral form DQ is given by 

DQ (27r)-P« dvolp(A) dvolg(S) drd^ Det«(^ - gB-'^t) V)e\P{B - tA-^g) , (7) 

where (ivolp(A) is a Haar measure for U(p) and dvolq{B) is the correspondingly normalized invariant measure for 
Herm^((7). Invariance of dvolq{B) means invariance under the transformation B gBg^ for any g G GLq(C); the 
explicit expression for dvo\q{B) is 

dB ' 

civol,(i?) ^ , di3=[]rfi?^^ n 2d9le(i?'-')dMS'^''') • 

^ ^ \=1 l<i^<iy'<q 

The symbol drda is short for the product of partial derivatives 



drd. = n n 



We mention in passing that the Berezin form ([7]) is invariant under the action of a Lie superalgebra 5i-{p\q) x Ql{p\q) 
where the first factor acts on the left of Q and the second one on the right. This invariance property actually 
determines DQ uniquely up to multiplication by a constant. The supermanifold of the supermatrix Q belongs to one 
of the ten families of Ricmannian symmetric superspaces described in2^. We also note that in the important case of 
p = q the Berezin form ([7]) is fiat, i.e., is given by a product of differentials resp. partial derivatives: 

DQ = {2n)-p' \dA\ dB d^d^ . (8) 

The superdeterminant is the usual one, satisfying lnSDct((5) = STr InQ with STrQ = Tr_B — TrA (bosons count as 
plus, fermions as minus): 

SDet"(Q) - (9) 

We repeat once more that the inequality n> q has to be satisfied in order for the formula ([6]) to be true. 

For the case that / is a Schwartz function (i.e., decreases faster than any power), a mathematical proof of the 
superbosonization formula ^ has been given in^^. More precisely, the theorem stated and proved in2^ assumes 
(besides n > q) that / is a holomorphic U(n)-invariant function of the vectors S'i , 5^ , Xi ; Xi (* = Ij ■ • ■ )'^)- These 
assumptions are always satisfied for the case of Gaussian disorder distributions. 

In summary, compared to previous works which suggest the existence a formula such as ([6]) , the main advance here 
is that the integration domain and sufficient conditions of validity have been specified. In the light of the discussion in 
the previous subsection, let us emphasize once more that there exist many choices of function F (although in practical 
applications there will usually be a natural choice) so that F{Q) becomes f{J2 V'i ^i) on substituting Q ^ . 
The theorem in^^ states that the superbosonization formula ^ holds true for any such choice of F and that paper 
also gives analogous formulas for the cases of orthogonal and symplectic symmetry. 

D. Proof of the superbosonization formula for p = q = n 

We now elaborate on the case of p = g = n, which will be a good starting point for making the generalization to 
n < q in Section III El An application of the special case ofp = q = n= l will be given in Section IIIII 

We simplify the notation by putting r := p = q = n. Thus we are now dealing with r vectors Si, . . . , Sr each of 
which has r complex components, and with r vectors Xit ■ ■ ■ iXr having r anti-commuting components. Let / be some 
function that depends only on the U(n)-invariant combinations (n — r) 

r r r r 

i—1 i—1 i=l i—1 

viewed as the blocks of a supermatrix. We wish to compute the integral of such a function: 



/.(2.)-^ / (10) 
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For this purpose we choose some function F{Q) of a supermatrix Q so that making the substitution 

(gill gsit). 

we recover from F the given function /. In the following we will change the integral / / in (jlOp to an integral over F 
in four steps. 

1. First step 

We begin by assuming that the vectors Si, . . . ,Sr are linearly independent, so that the rxr matrix which is formed 
by taking the vectors 5*1 , . . . , S'r to be the columns of that matrix is regular or of full rank. (The condition of linear 
independence is not always satisfied, of course. However, the sets of linearly dependent vectors form a set of measure 
zero. They can therefore be ignored for the purpose of integration.) Then we define square matrices a and r of size 
rxr and with anti-commuting entries by 

r r 
i=l 1=1 

By the chain rule of differentiation we have 

d^d^ = Dcf S*, 5,) drda , 

where drda- = Y\d^ / dr^^^dcr^'^ . By the assumption that the vectors 5*1, . . .,Sr are linearly independent, the rxr 
matrix 

r 

B■.= J2S^<»S^ 

i=l 

has non-vanishing determinant and hence an inverse. Thus we can write 

i=l 

and our integral (jTU]) then becomes 

/ = (27r)-''' J dSdSDef'iB) drd„ F (^^^ ^ 

2. Second step 

Now we take the Hermitian (and non-negative) matrix B — Si iS> Si as the new set of commuting variables of 
integration. Let d{B) be the (5-distribution (or (5-function) centered at zero on the linear space of Hermitian rxr 
matrices B, and let 6{B — J^^i ® ^i) shifted ^-distribution. We then have 

1 = J dB6(^B-J2St(^S^^ , 

where the integral is over the linear space of Hermitian rxr matrices B with flat (or translation-invariant) measure 
dB = Ylk dBkk Y[i<j 2 d'i)\z{Bij) d3m{Bij). Now if /i is an (integrable) function of ^ 5^ (X) 5^ , then by inserting the 
relation above and changing the order of integration we get 

J dSdSfi S*, (g) S*,) = J dSdS (^J dBS (^B -J2S^<E) 5,)^ A S^ (g) S, 

= J dBfi{B) J dSdSS(^B - ^ 5, ® S, 
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To compute the integral J{B) := J dSdSS (-B — (g) Si), let g G GLr(C) be any invertible r x r matrix with 
complex entries, and consider 



J{gBg^) = J dSdS S (gBg^ _ ^ 5^ ® 5^) . 



Of course, since '^Si® Si > , the function J{B) vanishes if one or several eigenvalues of B are negative. Changing 
variables from Si to gSi {i = 1, . . . , r) gives 

JigBg'^) = j dSdS Deflff^t) 5(g{B-Y,Si® S^) g^) . 

Now the (5-function obeys the transformation rule 5{gBg'^) = Dci~^ [gg'^) 5{B). We therefore obtain J{gBg^) = J{B). 
Since the action of GLr(C) on the positive Hermitian matrices hy B ^-^ gBg^ is transitive (i.e., every positive Hermitian 
matrix B can be written in the form B = gg^), it follows that J{B) is a constant independent of B > . 

This constant is readily found by taking B to be the r x r unit matrix, B = \r ■ Contributions to the integral 



J{lr) = j dSdS 6 [ir -^S^®S, 



come only from those S\,...,Sr that form an orthonormal system. Therefore J(lr-) = volU(r). More precisely, if we 
make the factorization {S\,S2-, . . . ,Sr) = e^e^^ with Hermitian X and Y, then dSdS becomes dB dk where B = e^-^ 
and dk for k = e'^ is the Haar measure for U(r) normalized by 

dk\Y=o = = n n 2 d^itiYij) ddm{Yij) . 

Thus J(lr) = J dk J dB 5{lr — B) = volU(r). After some of calculations (passing through the eigenvalue representa- 
tion of A:) one can show that with this choice of normalization the total volume of U(r) is 



volU(r) = / dk 



I 



(27r)K'-+i)/2 

U(r) 



0!1! •••(r - 1)! 

U(r) 

Altogether we obtain 

j dSdS fi (J2 S^ 5,) = volU(r) j dB fi{B) 

B>0 

Applying this relation to our situation gives 



/ 



/= (27r)-'''volU(r) J dBBet''{B)drdaF (^^^^^ '^^ 

B>0 



3. Third step 

We now re-express the volume factor volU(r) as a dummy integral volU(r) = /u^^j dk . Having done so, we view 

U(r) as a real subgroup of the complex group Ghr{C), and we extend the Haar measure dk (regarded as a differential 
form) from U(r) to GLr(C). For example, in the case of r = 1 this means that we set k = e'^ with y G [0,27r] and 
extend the given Haar measure dy of the unit circle U(l) C C to the holomorphic differential form dy = {ia)~^da on 
GLi(C) by letting a = e'^ be any non-zero complex number. 

In the general case, if Aij are the matrix elements of ^4 e GLr(C), the Haar measure dk extends from U(r) to 
GLr(C) as the holomorphic differential form 

where the variables Aij for i, j = 1, . . . ,r are regarded as a set of independent complex coordinates subject only 
to the condition Det(y4.) ^ 0. (Mathematically speaking, by interpreting the Haar measure as a differential form we 
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choose a fixed orientation of the unitary group.) To verify this formula for dk, notice that Det~^{A)dA is invariant 
under the transformation A i— > gAh for any g,h €z GLr(C). This invariance property in fact determines the Haar 
measure dk uniquely up to multiplication by a constant. To verify the normalization constant, one sets A — e'^ and 
notices that Det"''(A) dA\Y=o = i^'^dY. 

Let us now assume that F ^ ^ j has been chosen as an analytic function of the left upper block. Then, 



using the fact that 



\ T B 



U(r) 



dk kij^j-^ki2j2 ' ' ' — 



holds for any Z > 0, we can write 

Jf ^ (27r)-'-' J dBDet'iB) J dk drd^F ^ ""^'^^ 

B>0 U(r) 

= (27ri)-''' j dBDe^iB) J dADet-''{A) drd^F 



B>0 U(r) 



4- Fourth step 

In the final step, we use the fact that the holomorphic form dA is invariant under translations. We exploit this 
invariance to make the (nilpotent) translation A ^ A — ctB^^t. Since the unitary group is a closed manifold (i.e., 
has no boundary), such a translation does not give rise to boundary terms. Therefore, the result after translation is 

S>0 U(r) 

To bring this result into standard form, we write 

^« = ^'-^^'^ (d^) (pi#(I)) ^^^^Det^(A)Dct^(5) , 

where the second factor (returning from differential forms to measures) agrees with our Haar measure dk for U(r), 
and the first factor is a correspondingly normalized invariant measure for Herm^ (r) . These were denoted in the more 
general situation of Section Hi CI bv d'vo\p{A) and dvo\q{B) respectively. Note that from the relation 

Det«(A - aB-^T) Det^iB - rA'^a) = Det''(A) Det''(B) 

for p = q = r, the present expression for DQ in fact coincides with ([7]). 

Since the ratio of determinants Det(-B)/Det(A — oB^^t) equals the superdeterminant of the supermatrix Q = 
A a' 



T B 



we now arrive at the desired formula 



I f = I DQSBct'-{Q)F{Q) , 



where it is understood that we integrate over the unitary matrices A £ \J{r) and the positive Hermitian r x r matrices 
B e Herm^(r). This completes the proof of the superbosonization formula ^ for p = q = n = r . 



E. The case of n < g 

We now turn to the case of p fermionic replicas and q> n bosonic ones. This case is not covered by the results o:^ 
and needs separate treatment. As before, we wish to calculate the integral (O. 

Let us first verify by inspection that formula ^ cannot be true in the present situation. For this purpose, let 
p — q — r (but r > n) for simplicity, and consider some function / that depends only on the combinations Xi®Xi 
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and X]r=i Si Si . After superbosonization the integrand F is only a function of A and B (and does not depend on 
a and r). As we noted above, the invariant Berezin form ([7]) for p = q — r is 

DQ = {2tt)-''^ dYO\r{A) dYOlr{B) drda'Def {A)Bct'' (B) . 
Also, using the expression Q for the superdeterminant, our integrand becomes 

SDet" (g)F(Q) = (27r)-'^'dvol^(A) Det'-"(A) dvol^(B) Det''+"(B) 9^9^ Det" (l^ - B'^rA-^a) . 
Now the Fermi integral 

drd„ Det" (1^ - B-^rA-^a) 

vanishes identically, since Det"'(lr — B~^tA^^(t) can be at most of degree n x 2r in the matrix elements of a and 
r, whereas the differential operator drda is homogeneous of higher degree 2r^ > 2nr . (Thus there are not enough 
anti-commuting variables in the integrand to give a non-zero result when taking all partial derivatives.) 

However, if we scale A and B out of Det"(lr — B^^rA^^a) by sending, say, a Aa and r — > i?r, then the 
-B-dependence of the integrand becomes 

dvol^(B)Det"(S)i^(-,S) ^dBDct''-''{B)F{-,B) , 

which gives rise to a singularity when one (or several) of the eigenvalues of the positive matrix B approach zero. This 
singularity is non-integrable if the integrand F goes to a non-zero constant in the same limit. Thus the right-hand side 
of the superbosonization formula ([6]) is ill-defined (of the type of x oo) in the present case. Based on the treatment 
given for n = r in Section III D( we will now derive the correct formula for n < r . 



1. First step 



We now have n vectors Si, . . . , Sn, each with q complex components. We arrange these as a q x n rectangular 
matrix S :~ {Si, . . . , Sn) so that 

n 

SS^ ^Y,S^^S^. 

i=l 

We decompose the rectangular matrix S as 




where a is an n x n square matrix while the block b is rectangular of size {q — n) x n . For generic S the square matrix 
a is invertible and hence is an element a = g of the group GL„(C). Defining Z ba^^ we then have 

By simple power counting one sees that the volume element transforms as 

dSdS = dZdZ dgdg'Det'^"\gg^) . 

Here our notational conventions are the same as before, i.e., dgdg := Y[ij=i '^d,d\c{gij) d'3m{gij) and the same expres- 
sion goes for dZdZ. We now make a further change of variables 5^(1 + Z^ Z)^^^^h , which results in 

To explain the above factorization of S let Mat^ „ (C) denote all complex q x n matrices with full rank n . Every 
S E Mat^ n(*^) '^^^ be decomposed as S* = Th with h G GL„(C) and T e Mat'^ n(C) being the truncation of a unitary 
qx q matrix to the first n columns (i.e., the last q — n columns are deleted). This decomposition is not unique. Indeed, 
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if = Th is such a decomposition, then S = {Tk ^){kh) is also such a decomposition for any k £ U(n). Thus the 
correct mathematical statement of factorization is 

Mat;„(C) = (U(g)/U(g- n)) Xu(„) GL„(C) . 
Locally - more precisely speaking: whenever Det^'^(l + Z"^ Z) 7^ - we can make a definite choice of k by taking 

^ - (ia +^fi)-V.) ^ Mat;,„(C) . 

The factorization S = Th then means that we regard S as being given by its regular part h e GL„(C) and a point of 
the Grassmann manifold lJ{q)/\J{q — n) x U(n), which is the set of realizations of C" as a unitary subspace of C^. 
From the standard fact that the metric tensor of l]{q)/\J{n) x U(g — n) can be expressed by 

Trdz(i„ + z^zy^dz^ {lg.„ + zz^)-^ , 

one easily finds that dZdZ Bet~'^ {1 + Z'^ Z) expresses the invariant measure of U(q)/U(q — n) x U(n). We henceforth 
denotes this invariant measure by dvol(T). Thus we can summarize the discussion of this subsection by 

S = Th, dSdS = dvol(r) dhdhBef^^^ikh^) , 

which results in the formula 

J f = {2TT)~P" J dvol(r) J dhdhDet''-"-{hh'i)d^dxF 

Here x is the p x n rectangular matrix x = (xi? ■ ■ ■ iXn), and x is the corresponding n x p rectangular matrix whose 
rows are the row vectors xi, ■ ■ ■ ,Xn ■ In other words, we have J2Xi ^ Xi — XX I'YliXi® Si = X'S'^j etc. 



( XX xh^T^ \ 

\Thx Thh^ry ■ 



2. Second step 

In the next step we take the matrix elements of a := xh'' and r := hx as our new anti-commuting variables. This 
gives (cf. Section Hid ip 

J f = {27:)-P^ J dvol(r) J dhdhBeiP+'i-^{hh^)drd,F(^^^^^l'^^ tWt^ " 

Then we make a change of (commuting) variables to B hh'^ . By the same reasoning given in detail in Section 
HID 21 we obtain 

j / = (2^)-P"volU(n) j dvol(r) j dBYieiP+''-^{B)drd.F(^^~^^ . 



B>0 



The following steps are also similar to before. We introduce a dummy integral over U(p) with (volU(p)) ^ x 
/u(^)Det-P(iA)dA^l: 

I ^ - ^^olW) I '"""'^^^ J J Det^(iA) ^^^^ ^ [ Tr TBT^ ) ' 

S>0 U(p) 

Then we make the shift A ^ A ~ aB^^r. Expressing the result in terms of the invariant measures dvolp{A) = 
Det"^(iA) dA and dvol„(B) = Det"" (B) dB , we obtain 

B>Q U(p) 

This formula is arguably more complicated than The reason for this is that the rank of ^27=1 '^i ® never 
exceeds n , so there exist at least q — n zero eigenvalues and the range of ^ S'i (8) Si is a submanifold of the boundary 
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of Herm^((j'). Our derivation parameterizes this submanifold by two factors. The first factor, T, describes the n- 
dimensional complement of the kernel space of ^ S'i (E" Si in C ; thus it describes the complex n-plane spanned by the 
vectors Si, . . . ,Sn in C^. (The set of such subspaces C" ^ C is in one-to-one correspondence with the symmetric 
space \J{q)/\J{n) x \J{q — n)). The second factor, B, is the operator J2Si(^ Si restricted to its complex n-plane C". 

The merit of the result (|lip is it that allows us to make an exact transformation of the original problem to 
supermatrix variables in the parameter range n < q . Unfortunately, the invariance properties of the integral on the 
right-hand side are not very transparent. In other words, the Lie superalgebra Qi{p\q) x acts as first-order 

differential operators on the functions / and F, and while the transformation behavior with respect to this action is 
very clear on the left-hand side, it is not easy to see how the desired behavior emerges on the right-hand side. 

F. Supermatrix model 

Let us now return to the formulation of Section III C[ which is valid for n > q . The superbosonization formula ([6]) 
allows us to replace the initial problem of computing the correlation functions for a Gaussian ensemble of Hermitian 
random matrices with variance matrix Cy , Eq. with the problem of computing the generating function 

» M 

Z[J] = HdQ, eM~F[Q]) > F = Fo + Fi , 
i=i 

M M M 

Fo [Q] = hJ2 '^^3 sQ.s + nY, STr In - iF ^ STr Q,s , (12) 

M 

Fi[Q] = -i^ STrQ, s (£ - F • 1 + J,) , 
j=i 

where E is the "center of mass" of the energy parameters in £. This is the most general supermatrix model for 
ensembles of Gaussian random matrices with unitary symmetry U(n)^^ {M — N/n). Note that no approximations 
have been made, and (|12p is an exact reformulation of the original problem. The entries of the variance matrix Cij are 
required to be real symmetric and positive, but are otherwise arbitrary. The size of the supermatrices Q depends on 
the correlation function to be calculated and is the same as for the supermatrices Q of the standard non-linear sigma 
modelii^. However, while the eigenvalues of the usual sigma model field Q are constrained to be ±1, the eigenvalues 
of our superbosonization field Qi fluctuate; they are real and positive in the boson-boson sector and unitary (i.e., of 
unit modulus) in the fermion-fermion sector. 

Let us also point out that, since the key formula (O holds for a large class of functions / (not just the Gaussian 
functions), the present method is not restricted to Gaussian disorder distributions. In the case of a more general 

disorder distribution, the term exp (^~^ X]fj=i Qi ^ Qj^^ p2p is replaced by a functional which is determined 

by the Fourier transform of the disorder distribution function. 

It should also be clear that such a description as p2| exists even for n < q . We just need to replace the integration 
domain and integration measure DQi SDet"((3i) by the modified one constructed in Section fll El 

III. DENSITY OF STATES FOR ALMOST DIAGONAL MATRICES 

We now demonstrate how the method developed in the previous section works for the density of states of almost 
diagonal random Hermitian matrices. We focus on the n = 1 case, which was investigated in^^, and we will compare 
the results of this reference with ours. Note that the limit of almost diagonal random matrices is not accessible (for 
n = 1) via the standard non-linear sigma model for random matrix problems. 

To calculate the density of states from the generating function p2p . we set p = q ^ 1 and s = 1. Our super- 
bosonization field Qi is now a 2 x 2 supermatrix 




with Berezin integral form 



DQi — {2ni) ^ dbi dui d^dai 
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The density of states per unit length for a system of N sites is expressed as 

N 

p{E) = (27riV)-i 3m ^ i(Tr Q,) , (13) 



i=l 



where (. . .) means the average with respect to the statistical weight e"^!*^! of p2|) . taken with vanishing source term 
Ji = . The diagonal matrix of energy parameters is £ = diag(£', E) with 3m E' > . 

We now proceed by first solving the diagonal variance matrix = Co Sij case exactly and then expanding in 
the off-diagonal terms C^ , i ^ j. For notational simplicity we assume that the system is one-dimensional with 
translation-invariant C^ = c(|i — j|) (although neither assumption is really necessary). 

In the zcroth order expansion the integral ^(TrQ^) = (6^) = (oi) in (|13p factors as a product of iV independent 
integrals, one for each site. The — 1 integrals for the sites j ^ i all are unity: 

/.0,SDe.«3,),— = 1 (C, = C,). 

which is a consequence of supersymmetry. The remaining integral for the distinguished site i is 

y i?Q,SDet(Q,)a^e'^STrQ.-(Co/2)STrQ,Q, 

= (27ri)-i / dbi daa,a,(6-Ta-V)e'^(''-'^)-«^«/2)('''-'-2--) = / e'^'''(^°/2)6^d6 , (14) 
Jm.+ "'U(i) ^0 

where in the second line we dropped the index i from the integration variables. The err term in the exponent cannot 
contribute to the Fermi integral drd^ , as one must pick the term —Ta~^a in front of the exponential in order to have 
a non-zero integral over a G U(l). Hence, denoting the density of states in the zeroth order of this expansion by 
Pd{E), we have (for 3m E -f Q+) 

PdiE) = n-'3mi{a,) = (27r)-i / e''"'~^^'>/2)b^ = (2^Co)-i/2e-^'/2Co (^g) 



This Gaussian function is of course none other than the probability distribution function of the diagonal elements of 
our almost diagonal random matrix. 

Now we calculate the correction Spd{E) coming from the off-diagonal elements. For an almost diagonal variance 
matrix Cy we may approximate the generating function Z[J] by expanding the exponential in the off-diagonal terms: 



N 



Z[J] = [llDQ. SBetm eiSTrQ.(£+./.)-(Co/2) STr q.q. (i _ 1 ^ STr Q,g, + . 

Therefore the correction to Pd{E) may be written as 

6p{E) = -(4^iV)-ijmi^^C,fe (Tr(Q,)STr(Q,Qfc)) . (16) 

Again, (. . .) is a product of N independent integrals, — 3 of which are unity by supersymmetry. Using the relation 
(Qf^^) = i(Tr Q) S^" and STr 1 = , we see that (. . .) in p6)) vanishes unless i = j oy i = k . Hi = k ^ j then (note: 
2 angled brackets included in r.h.s.) 

(Tr(QO STr(Q,Q,)> - i(Tr(QO STr(Q,)) (Tr(Q,)) . 

The single-site average of TrQ^ is given by (fT4|) . The expression for the single-site average of Tr((5i) STr((3i) is the 
same except that an additional factor STrQ,; — bi — ai has to be inserted under the integral sign. The term — of 
this factor gives vanishing residue at the simple pole = and hence does not contribute. So only an extra factor 
of bi = b remains and we get 



J i?Q, SDet(g,)Tr(g,)STr(QOe'^ STY Q.-(Co/2) STr Q.Q, ^£ 



Altogether we then obtain for the correction to Pd{E) (when the variance matrix C^ is almost diagonal): 

5p{E) = -^^ (e-^VCoerfi(£;/V2^)) , (17) 
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where Ci = X^j^i ^ij erfi(x) is the imaginary error function, 

Note that Sp{E) is a total derivative. This is as expected as turning on the off-diagonal elements of C,j neither creates 
nor destroys levels but just changes their positions. 

The high-energy limit of the density of states for almost diagonal random matrices was calculated in Ref.i^. To 
compare with those results, we note that the large- a; limit of the imaginary error function is x e~^ erfi(a;) — > ^/n , 
which gives 

SpiE) = ^Pd{E) [E » . 

Thus we have the simple result that the density of states for an almost diagonal Gaussian random matrix at high 
energies (i.e., in the far tail of the Gaussian distribution) is 

p{E) = ^(Q,/C,,) pd{E) {E » . (18) 

j 

This result is equivalent to that calculated in Ref.— . 



IV. REDUCTION TO THE STANDARD SIGMA MODEL 



Let us now see how the generalized supermatrix model (|12|) reduces to the standard diffusive non-linear sigma model 
under suitable conditions. In order for such a reduction to take place, the superbosonization field must get localized 
in a certain low-energy submanifold (corresponding to the standard sigma model field) of the total field space. As we 
shall see, the latter comes about if the number of orbitals n is large. If, in addition, the superbosonization field has 
enough collectivity (or stiffness) due to a variance matrix Cij with sufficiently long range, then the effective degrees 
of freedom of the problem are the Goldstone modes associated with the low-energy manifold, and one recovers the 
diffusive non-linear sigma model. 

To identify the low-energy manifold, we must first understand the symmetries of the functional -Fb[Q] in (I12p . 
Beginning with the boson-boson sector, let g € GLg(C) be any invertible complex q x q matrix and consider the 
transformation 

1-^ Ai , Bi gB^g^ , ai i-> aig^ , n gn . 
This transformation is a symmetry of the Bcrezin integral form DQi of ([7]) . It is a symmetry of Fq [Q] if 

where Sq is the boson-boson part of s . The condition (I19p singles out a pseudo- unitary subgroup of GL^ (C) . This 
non-compact group is denoted by U(g+, q^). 

Turning to the fermion-fermion sector, let (g, h) e U(p) x U(p) and consider the transformation 

Ai i-> gAih^^ , Bi^ Bi , ai ^ goi , n nh^^ . 

Again this is a symmetry of DQi . In order for it to be a symmetry of Fq [Q] we must impose the condition g = h , 
which singles out the diagonal subgroup U(p) C U(p) x U(p). 

The product of groups U(p) x U(g+,(7_) is the group of bosonic symmetries of our problem. There also exist a 
number of fermionic symmetries. It is not difficult to show that on inclusion of these symmetries the symmetry group 
becomes a Lie supergroup U(p | q^,q^). This is known to be the symmetry group of the standard non- linear sigma 
model (for systems in the unitary symmetry class). Note that the group action of U(p|(7+,(7_) is by conjugation 
Q,s^UQ^sU~\ 

The next step is to look for minima of the 'energy' functional Fo[Q] in p^ . (The term Fi is considered small for 
present purposes). Varying this functional gives the saddle-point equation 

= Jq-Fo[Q] - C^J^Qo^ - iQi)~^n - iEs . (20) 



7^505 = So 



So 
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Assuming translation invariance of (and setting C := Co + Ci = C^) one first looks for i-independent solutions 
Qi = Q m the space of diagonal matrices. From (pO]) each eigenvalue of the diagonal matrix Q satisfies a quadratic 
equation, which in general has two different solutions, 5'^+^ and q^~\ In the boson-boson sector, one of these solutions, 
say q'^~\ is ruled out by the positivity condition > . However, in the fermion-fermion sector both solutions g^^-* 
are in principle admissible. Connected supermanifolds of solutions of the saddle-point equation are then generated by 
the conjugation action Qs ^ UQs U^^ of the symmetry group. The saddle-point manifolds thus obtained are orbits of 
\J{p I g+, g_) which can be classified by the number of eigenvalues q^^^ and q^~^ of the solution in the fermion-fermion 
sector. For energies 'inside the band', \E\ < \/ 4Cn , all these supermanifolds are parameterized by 

Qs = C-^/^{n^ E^/4Cy/^UAU-^ +i{E/2C)l , U E V{p\q+,q-) , (21) 
where A is a diagonal matrix with eigenvalues ±1 : 

A = diag(si, So) ■ 

The boson-boson part of A is uniquely determined by the boson-boson part sq of s , but the signs in the fermion- 
fermion part si of A are arbitrary. In the general situation, each of these saddle-point manifolds (corresponding to 
different choices of Trsi) make a contribution to the generating function ()12|) . 

In the following we focus on the important case p — q . By making the maximally supersymmetric choice si = sq , one 
obtains a distinguished saddle-point manifold which dominates (under conditions to be specified) , while contributions 
from other saddle-point manifolds are suppressed by fermionic zero modes due to a breaking of supersymmetry. 
Low-energy configurations of the superbosonization field are now expressed as 

Q,s = TTvQ,+i{E/2C)l , V ^-K-^C-^/^in- E"^ /ACf^ , (22) 

where Qi is a dimensionless field, namely the standard sigma model field 

= [/^AC/-\ A = diag(so,so) , e U(p |p+,p_) . 

By computing the second variation of Fq [Q] at the minimum, one sees that the fluctuations of modes transverse to 
the low-energy manifold (the so-called massive modes) are controlled by the quadratic form 

hij = {-Kvfdj + nSij n{Sij + C^j/C) . 

(Here the second expression makes the simplifying assumption that E — .) Since this quadratic form is bounded from 
below by the diagonal form n 6ij , its eigenvalues are never smaller than n and it therefore follows that fluctuations of 
the massive modes are strongly suppressed in the limit of n ^ 1 . 

In the case of small n, when is short ranged (i.e., our random matrix is almost diagonal), the massive modes 
fluctuate strongly and there is no controlled reduction to the low-energy manifold of flelds Qi . One might now think 
that the situation gets better when Cij is taken to be long ranged. However, this is definitely not true when the 
superbosonization method is used! The problem is that the Fourier transform C(fc) of a long-ranged variance matrix 
Cij ~ c{\i — J I) decreases with increasing wave number k and is close to zero at high wave numbers (near the edge 
of the Brillouin zone), which implies that our high- momentum massive modes always have a small mass (of order n) 
when n is small. 

It should be stressed that the case of a long-ranged variance matrix is much better handled by the traditional 
Hubbard- Stratonovich transformation approach. There, fluctuations of the massive modes (say at i? = 0) are con- 
trolled by another quadratic form, 

h^f^^ = n{Sij + Cwij) , 

where Wij is the matrix inverse of Cij (i.e., '^j WijCji = Sn). The Fourier transforms of Cij and Wij are reciprocals 

of each other, iv{k) = C{k)~^. Thus, if the eigenvalues of Cij become small near the edge of the Brillouin zone, then 
those of Wij become large, and therefore the high-momentum massive modes of the Hubbard- Stratonovich approach 
are truly massive for the case of long-ranged Cij (and all n including n = 1) and can be integrated out in a controlled 
way. In summary, since a different quadratic form enters the game, the reduction to the standard sigma model in the 
present approach is valid under conditions that are not identical to those of the Hubbard-Stratonovich approach. 
We now insert ([22]) into the expression (fT2|) for Fq to obtain an effective energy functional for the Qi fleld: 

i,3 
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This is a good approximation when n ^ 1 , in which case we may simply neglect the massive modes. On the other 
hand, if n is small the massive modes fluctuate strongly and we face the non-trivial task of integrating them out with 
a non-perturbative calculation. While this will not change the symmetries of the effective energy functional for the 
Qi field, it may cause a major renormalization of the coupling parameters Cij . 

Therefore, to maintain quantitative control of the situation, we now assume n to be large. Taking into account the 
term Fi in (fT2|) . we arrive at a low-energy effective action 

M 

F[Q] = 5 ^ STr Q,Qj - ini^ ^ STr Q, {£ + J,) , (23) 

i^j 1=1 

where diagonal terms have been dropped from the double sum over i,j because STr QiQi — . Recall that Qi is the 
standard sigma model field. 

The action (|23p describes a spectrum of physical situations ranging from strong localization to diffusive behavior. 
Strong localization occurs in the limit of almost diagonal random matrices, where our coupling coefficients Cij {i ^ j) 
are small, leading to a disordered sigma model field Qi . (Note that small Cij can be realized in spite of n being large.) 

The diffusive regime can be realized by taking the variance matrix Cij to be long ranged. In this limit, spatial 
variations of Qi are suppressed and we may make a continuum approximation Qi — > Q{r), expanding the first term 
on the right-hand side of in gradients of Q{r). A standard computation gives the diffusive action 

i STr Q,Q, « - ^ STr ^ dr D (VQir))^ , 

where the 'diffusion coefficient' is D — 2^1/ J^j N ~ J p^'y . To summarize, the discussion above has shown how, in 
some cases, the conventional non- linear sigma model can be obtained by reduction from the supermatrix model (jl2p . 

Now in various regimes such a reduction cannot be done and one has to work with (fT^ . Among these are models 
with a critical point where an Anderson metal-insulator transition takes place. Much attention has been given to such 
models^ii^i^^'^ with the aim of better understanding the critical behaviour at the transition, ft is an interesting 
open question whether the supermatrix model (|12p can be of use in the analysis of these critical models. 

V. CONCLUSION 

The main outcome of this paper is an exact mapping of a Gaussian random matrix problem to a supermatrix model. 
The measure of integration and the structure of the matrix field Qi are given in Eq. ([7|) for n> q and in Eq. (fTTj) 
for n < q . One might be tempted to call the supermatrix model a 'generalized non-linear sigma model', but such a 
terminology would be misleading as the target space of the supermatrix model is not homogeneous (unlike with what 
is called a non-linear sigma model in the strict sense of the word) with respect to its symmetry group. 

A notable result is that the superbosonization formulas ^ and (jlip are different, depending on the value of the 
number of 'orbitals' n relative to the number of commuting variables q in the field ip . The former determines the 
real space structure of the random matrix ensemble, while the latter is equal to the number of points q in the (/-point 
correlation function in the energy space. Thus the superbosonization formulas (O and pip indicate that there exists 
a correlation, in a certain sense, between real space and energy space. This fact seems to be fundamental, although 
its implications are not yet fully understood. 

In this context, let us make one observation which we consider to be of relevance. It has been a puzzle for many 
years now how an inspired use of the replica trick with fermionic replicas correctly reproduce a^^'^^i^^'^" the DoS-DoS 
correlation function in the large-n random matrix limit, whereas the replica trick with bosonic replicas is known to 
fail^ when the same limit is invoked. Based on the difference between q < n and q > n observed above, we propose 
the following resolution of this long-standing puzzle. 

In the replica trick one needs to calculate the observable of interest for every number of replicas and then, using 
all this information, one tries to find an analytic continuation to zero replica number. When using the replica trick 
with q bosonic replicas, we conjecture that vital information about the analytic continuation g is contained in 
the high replica numbers q > n , where the behavior of the correlation functions seems to be qualitatively different 
(based on what we have seen with the superbosonization formula) from that for q < n . From this vantage point, it 
has to be regarded as an ill-advised scheme to take the large-n limit (and thus the saddle-point approximation leading 
to the bosonic replica NLcrM) at an early stage of the calculation, as was done in^^. Indeed, in the process of taking 
n — !■ GO the q > n branch, and hence all information carried by it, is lost from the computation. Our explanation of 
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the long-standing puzzle is that it is the loss of the q > n branch which makes it impossible to construct the correct 
analytic continuation to zero replica number. 

This scenario and other questions related to the difference between the cases n > q and n < q are interesting 
directions for future research. Here, in contrast, wc outline the idea of a uniform representation of both cases which 
requires a non-trivial limiting procedure. For simplicity we restrict ourselves to the case otn = l, p = q = 2. As 
was explained in detail at the beginning of Sect. Ill El the simple superbosonization formula fG]) fails in this case. For 
the sake of discussion let us recall the details here. The argument is especially transparent when the function F{Q) 
depends only on the commuting variables in the matrices A and B and is non-zero in the limit B ^ 0. Note that the 
left-hand side of ([6|) is non-zero in general. On the other hand, since A, B, a, and r for p ~ q = 2 are 2x2 matrices, 
the superdeterminant 

SDet(g) = Dct(B) Dct"^(A) Dct(l - A" V) 

is a quartic polynomial in eight Grassmann variables cr'"^ and t"^^ [ji, — 1,2) and therefore Yl dr^^i^dai^'' SDet{Q) — . 
This means that the right-hand side of ([6]) is zero, unless there is a divergence in the integration over commuting 
variables. Here one must pay attention to the fact that the integral over the bosonic variables B can be singular 
along the Det(i3) = boundary. A quick way to see this is to make a transformation to new variables a' — A~^a 
and r' = B^^t , which has Jacobian Det~'^{A)Det~''{B) and leads to a singular dependence Det^~^(i3) = Det~^(_B) 
at q = 2 . Thus, a naive application of the superbosonization formula ([6]) to the case n^l, p = q = 2 leads to an 
ill-defined expression of the type x oo . 

The idea of a unified description stems from the observation that the deficiency of the Grassmann variables in 
SDet((5) (where 4 are missing from the full set of 8), which is responsible for the part in x oo, does not persist 
when SDet((5) is replaced by SDet^^"(Q) with a G C being a small regularization parameter. Of course, to move 
a off zero, we must first define the factor Det~^~°'{A) in SDet^'*'"(Q). This can be done by making some choice 
of fundamental domain for the logarithm In : GL2(C) C. Restricting A E U(2) C GL2(C) to this fundamental 
domain, we may expand 

SDeti+"(g) = SDet(Q) (1 -I- a InSDet(Q) ...) . 

One then easily verifies that the Fermi integral of SDet^^"((5) is non-zero for small but non-zero a. At the same 
time, the singular factor Det~^{B) in the integrand is replaced by Det"~^(i3), which makes the singularity at _B — > 
integrable for IHe a > . Moreover, it is possible to analytically continue the integral over A G U(2) to a 7^ by 
taking the integration contour, i.e., the 4-dimensional real submanifold of integration in GL2(C), to infinity along the 
boundary of the fundamental domain of the logarithm. (This will work if the boundary of the fundamental domain has 
been chosen as a submanifold of rapid decay of the integrand.) Thus, carefully taking the limit limc(^o+ SDet^^"((5) 
instead of immediately setting a = , is a way to give meaning to the ill-defined expression x 00 . 

Exploratory calculations done along these lines give a remarkable coincidence with Eq. (|lip . Our explorations 
suggest that the superbosonization formula in both the n — q — 1 and n = 1, q ^ 2 cases can be written in a uniform 
way resembling the original form of Ref.— and Eq. ^ but with limQ^o+ SDct^^"(Q) replacing SDet{Q). It should 
be stressed, however, that the mathematics at hand gets more complicated as the number of replicas p — q goes up. 
We have not yet made a serious effort to confront these complications and more work is needed to put the idea of 
a-regularization on a solid mathematical basis in the general case. 

To finish this Conclusion, we would like to compare our superbosonization (SB) approach with the standard one^i^ 
based on the Hubbard-Stratonovich (HS) transformation and the degenerate saddle-point approximation. As was 
mentioned in the Introduction, the main disadvantage of the standard approach is that it does not apply to random 
matrix ensembles with a short-ranged variance matrix Cij . Another disadvantage is that it is applicable only to 
Gaussian ensembles. Both restrictions are lifted in the new SB approach. The new approach is in some sense dual 
to the standard HS approach as it relies on the variance matrix Cij , while the HS approach relies on its inverse, 
{C^^)ij . This feature makes the new SB approach most efficient for short-ranged Cij , e.g., for almost diagonal 
random matricesi^. It can also be used for long-ranged Cij but only under the condition that the number of orbitals 
n is large. For n = 1, extracting the standard NLcrM from Eqs. ([12]) is a highly non-trivial task for the case of a 
long-ranged matrix Cij (which is a textbook example of the derivation of NLcrM in the framework of the HS approach). 

This example shows once again that the two exact supermatrix representations of the random matrix ensemble 
which emerge after the SB and HS transformations, are fundamentally different and largely complementary. Thus the 
new SB approach does not negate or supersede the standard HS one. Rather, it extends the supersymmetry method 
to an area previously not accessible to it. 
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